library(sasld)


# 9.1 ---------------------------------------------------------------
binom.test(x=10,n=10,p=0.5,alternative="greater")

binom.test(x=22,n=54,p=1/7,alternative="greater")
sum(dbinom(22:54,54,1/7))
# -------------------------------------------------------------------




# 9.2 ---------------------------------------------------------------
data(anorexia)

ok <- which(anrx$therapy == "CB")
y <- anrx$after - anrx$before; y <- y[ok]
t.test(y,mu=0,alternative="two.sided")$p.value

y[15] <- 2.9
t.test(y,mu=0,alternative="two.sided")$p.value

y[15] <- 20.9
for (i in 1:29) print(c(y[i],t.test(y[-i],mu=0,alternative="two.sided")$p.value))
# -------------------------------------------------------------------




# 9.3 ---------------------------------------------------------------
set.seed(123456)
x <- simProp(nsim=10000,n=116,p0=1/3,alternative="greater")
summary(x)
summary(x,alpha=0.01)


set.seed(654321)
simProp(nsim=1,n=116,p0=1/3,p1=0.5,alternative="greater")

set.seed(654321)
x <- simProp(nsim=10000,n=116,p0=1/3,p1=0.5,alternative="greater")
summary(x)
# -------------------------------------------------------------------